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ABSTRACT  ;  'The  axisymmetrlcal  elast 
shell  under  longitudinal  impact  is 

;ic  deformation  of  a  circular  cylindrical 

1  investigated  by  using  a  system  of 

nonlinear  equations  with  the  propagation  of  elastic  stress  waves  taken 
into  account,  and  without  any  assumptions  concerning  the  mode  of  buckling. 
One  end  of  the  shell  is  fixed,  the  other  end  is  axially  impacted  by  a 
rigid  solid  moving  at  a  velocity  Vj  the  ratio  m  of  the  mass  of  the  body 
to  the  mass  of  the  shell  is  given.  The  analysis  of  the  impact-deformation 
process  in  this  shell  is  reduced  to  solving  this  nonlinear  system  with 
initial  and  boundary  conditions  by  the  method  of  finite  differances, 
utilizing  an  explicit  scheme  whose  convergence  and  stability  was  checked. 
The  behavior  of  the  shell  was  studied  in  the  time  interval  in  which  the 
longitudinal  compression  wave  propagates  along  the  whole  length  of  the 
shell,  and  the  first  reflected  wave  comes  back.  The  results  from  cal¬ 
culating  the  normal  displacements  along  the  shell  at  various  instants  of 
both  waves  propagating,  for  the  ratio  m  =  5.64  and  nondimens ional 
velocities  V/a  =  0.0005;  0.001;  0.002;  and  0.004  (where  a  is  the  velocity 
of  sound)  are  shown  in  diagrams  and  are  examined.  The  qualitative  aspect 
of  the  shell  deformation,  especially  the  formation  of  maximum  local  dis¬ 
placements  during  the  passage  of  both  the  compression  and  the  reflected 
waves  as  related  to  V  is  discussed  and  found  to  be  in  agreement  with  the 
A.  Koppa  phenomenological  theory  based  on  experimental  results.  Orig. 
art.  has;  5  figures  and  15  formulas.  English  translation;  9  pages. 


TRANSCRITICAL  DEFORMATION  OF  A  CYLINDRICAL  SHELL  ON  IMPACT 

V.I.  Borisenko,  A. I.  Klokova 
(Kiev) 

On  the  baeie  of  nonlinear  equations  whioh  take  account  of 
shear  and  rotational  inertia^  the  axially  symmetric  elastic  de¬ 
formation  of  a  circular  cylindrical  shell  in  longitudinal  impact 
is  considered. 

The  deformation  of  the  cylindrical  shell  in  longitudinal  im¬ 
pact  is  investigated  on  the  baeie  of  wave  equations  without  re¬ 
sorting  to  complementary  assumptions  concerning  the  buckling  con¬ 
figuration  of  the  shell.  It  is  assumed  that  one  edge  of  the  shell 
is  supported  while  the  other  is  subject  to  longitudinal  axially 
symmetric  impact  by  a  rigid  body  moving  at  the  velocity  Fq. 

Experimental  studies  of  the  process  of  loss  of  stability  of 
a  cylindrical  shell  on  Impact  [5]  showed  that  the  wave  nature  of 
stress  propagation  in  the  shell  markedly  affects  the  nature  of  the 
loss  of  stability.  In  this  connection,  most  of  the  studies  of  the¬ 
ory  of  this  problem  [1,  3]  analyze  the  buckling  process  on  the 
basis  of  systems  of  equations  of  parabolic  type  which  make  no  al¬ 
lowance  for  the  process  of  wave  propagation  in  the  shell.  More¬ 
over,  the  solution  of  the  equations  normally  is  associated  with 
specified  buckling  configurations. 

To  investigate  the  process  of  the  deformation  of  the  shell 
on  impact  we  will  employ  the  equations  derived  by  M.P.  Galln  [4]. 
For  elastic  deformation  these  equations  are  given  in  the  monograph 
[7]  on  adding  terms  that  take  account  of  the  geometric  nonlinear¬ 
ity.  The  equations  pertaining  to  the  axially  symmetric  case  are; 


dx 

dx  R  , 


dM,, 


dx 


.Q, _ 

12 


where  the  notation  is  the  same  as  in  [7], 


The  stress-strain  ratios  are  as  follows: 


(1) 
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1  — v>[/? '^''[Sje'*' 2  ^ajcj  Jj*  12(1— v»)a;e* 

Equations  (1)  make  It  possible  to  Investigate  the  propaga¬ 
tion  of  longitudinal  and  transverse  perturbations  In  the  shell, 
which  cannot  be  done  by  proceeding  from  the  classical  equations  of 
the  shell  theory  based  on  the  Klrchhoff-Love  hypothesis.  These 
relations  were  derived  [*»]  from  general  nonlinear  equations  of  the 
theory  of  elasticity  by  expanding  the  displacements  Into  a  series 
with  respect  to  the  power  of  m  (the  *-axls  Is  directed  at  right 
angles  to  the  median  surface  of  the  shell)  and  preserving  two 
terms  of  the  series  for  tangential  displacements  and  one  term  for 
nozmial  displacement.  There  also  exist  other,  more  rigorous  methods 
of  reducing  to  two  dimensions  the  three-dimensional  problems  of 
the  theory  of  elasticity  [5]. 

In  the  formula  for  stress  Qi  from  (2)  the  coefficient  » 

■  0.86  Is  Introduced  In  accordance  with  [8]. 

Following  the  substitution  of  Relations  (2)  In  (1)  we  have 
the  system  of  differential  equations 

^  I  ,  ^fd^u  ,  V  dw\  ~  .dip  ..  i5«. 

d*9  ,  6  f  dw  _  1 

dt»  "  aa»  e»  \3o  ’ 

where  a  ■  */F,  x  ■  at/R,  3  *  w/R  are  a  dimensionless  coordinate, 
time  and  elastic  displacement,  respectively; 


u  =  ulR-, 


.  AMl-v) 

6 - 2 - • 


Assume  that  the  edges  of  one  side  of  the  shell  are  rigidly 
fastened.  Then  the  fastened  side  must  satisfy  the  conditions: 

u  a  w  SB  9  =  0,  (a  a  0).  ( ^  ) 

The  side  subjected  to  the  Impact,  on  the  other  hand,  must  satisfy 
not  only  the  two  conditions 

w  BB  <p  CB.  0.  (a  =  n),  ( 5 ) 

but  also  the  condition  that  will  obtain  If  the  force  of  Inertia 
acting  on  the  striking  body  Is  equal  to  the  longitudinal  stress 
acting  on  the  edge  of  the  shell 
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(6) 


M^  =  -2hRT, 

In  Eqs.  (5)  and  (6)  y  »  l/B  Is  the  dimensionless  length  of 
the  shell  and  M  is  the  mass  of  the  striking  body. 

Using  Relation  (2)  for  Tn  and  passing  on  to  dimensionless 
coordinates  and  displacements,  on  taking  account  of  (5)  we  have 

[!+?©■]• 

where  m  is  the  ratio  of  the  mass  of  the  striking  body  to  the  mass 
of  the  shell. 

Considering  that  prior  to  the  impact  the  shell  was  in  a  state 
of  rest,  we  write  the  initial  conditions  as 


^■=0;  0<a<y.  (t  =  0); 


du 

dx 


a  =  y,  (T  =  0); 


0  <  a  <  i. 


(T  =  0); 


(8) 

(T  =  0). 


The  investigation  of  the  process  of  deformation  of  a  rigidly 
attached  cylindrical  shell  subjected  to  a  longitudinal  impact  by 
a  rigid  body  reduces  to  the  solution  of  Eqs.  (3)  with  initial  (8) 
and  boundary  (4),  (5),  (7)  conditions.  In  accordance  with  the 
Cauchy-Kovalevskaya  theorem  this  problem  has  a  solution  that  ap¬ 
pears  to  be  unique.  Such  a  circumstance  makes  it  possible  to  in¬ 
vestigate  the  process  of  deformation  of  a  shell  lacking  any  ini¬ 
tial  irregularities.  The  solution  will  be  sought  in  numerical 
fonn  with  the  aid  of  the  net  method.  We  utilize  the  following  dif¬ 
ference  relations: 

N’t+l.n  — 
dp  “  2A 

fQi 

ap’  *  A’ 

where  ip  is  the  differentiable  function,  B  is  an  independent  vari¬ 
able,  A  is  a  mesh  of  the  net  with  respect  to  the  variable  B,  and 
tpy  j  is  the  value  of  the  function  at  the  net  node  numbered  t,j. 

In  the  region  of  the  variables  t,  a  we  Isolate  individual 
points  —  nodes  of  the  net  with  a  uniform  mesh  spacing,  as  is  shown 
in  Pig.  1.  The  mesh  with  respect  to  the  variable  a  is  denoted  as 
t  and  the  mesh  with  respect  to  the  variable  t,  as  h.  The  subscript 
i,  denoting  the  number  of  the  vertical  line  (a  ■  const)  ranges  in 
value  from  zero  to  k  (a  ■  y).  The  subscript  j,  denoting  the  num¬ 
ber  of  the  horizontal  line  (x  ■  const)  ranges  in  value  from  0  to 
Jf,  where  R  depends  on  the  Interval  of  time  over  which  the  solution 
is  sought. 
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Pig.  1 


With  the  aid  of  Relations  (9)  we  rewrite  System  (3)  as 

Ui./+1  =.2^1  —  ^  j  Uij  —  +  U/-I./)  T 

-r  ^  (“'/+!./  —  “'/-ij)  — 

—  ^(9/+u~‘P/_i./)J  +  /  j^‘P/+u  'P/-I./ 

-- r<“'^+M-2«'/./+“'/-u)]--  (10) 

B'/J+l  •=  (2  —  /l*)  +  -p  {Wt+\J  —  2k»/j  +  X 

^[ft  +  ^  (“/+!./  —  «/-!./)  +  j  +  ^  (W/+1./  —  Wl-\j)  X 

^  [t  '  "**  X  —  ®<-i./)  j  ~ 

6ft’  ,  t  ft’v  , 

2/  2/ 

'P/./+I  2?  j  '^‘•1  ~W 

^  6ft”*  Y 

I  2l(?  “'/-i./^' 

Thus  the  known  values  of  the  functions  u,  u  and  \|»  in  the 
{j  —  1)^  and  jth  rows  can  be  used  to  determine  from  Formulas 
(10)  the  values  of  these  functions  in  the  (j  +  l)th  row.  Such 
difference  schemes  have  become  termed  explicit  [217  As  is  known, 
the  advantage  of  explicit  difference  schemes  over  their  implicit 
counterparts  lies  in  the  simplicity  of  the  computational  algo¬ 
rithm.  However,  explicit  difference  schemes  do  not  always  display 
computational  stability.  This  stability  largely  depends  on  the 
selected  ratio  between  the  net  meshes  and  their  dimensions. 

Hence,  to  select  the  ratios  between  the  dimensions  and 
meshes  of  the  net,  the  linear  parts  of  difference  equations  (10) 


-  it  - 


were  subjected  to  a  numerical  analysis  on  a  digital  computer  with 
the  aid  of  the  so-called  e-scheme  [2].  As  a  result.  It  turned  out 
that  the  linear  part  of  difference  scheme  (10)  Is  unsuitable  for 
computation,  since  the  roundlng-out  error  rapidly  Increases  with 
transition  from  one  row  of  the  mesh  to  another.  One  of  the  rea¬ 
sons  for  this  Increase  In  error  Is  the  presence  of  high  coeffici¬ 
ents  In  the  right-hand  parts  of  Relations  (10).  To  reduce  these 
coefficients  the  substitution 


*F/./ —  2  (11) 

was  carried  out  In  the  third  formula  of  (10). 

Subsequent  calculations  on  the  basis  of  the  e-scheme  showed 
that  In  the  presence  of  e  the  coefficients  increase  Insignificant¬ 
ly  with  transition  from  one  row  of  the  nodes  to  the  next.  Further 
verification  of  the  convergence  and  stability  of  the  difference 
scheme  consisted  In  that  the  necessary  calculations  were  carried 
out  on  successively  reducing  the  net  meshes  until  the  findings  ob¬ 
tained  before  and  after  the  reduction  of  the  mesh  coincided  to  ' 
specified  degree  of  accuracy. 

With  the  aid  of  initial  conditions  (8)  we  determine  the  val¬ 
ues  of  the  sought  functions  In  the  first  two  horizontal  rows  of 
the  net  (cf.  Pig.  1) 

j  =  0;  u,.  g  =  ,,=*  <P,  0  =  0,  (i  =  0,  1 , . . . ,  A); 

/  =  1;  =  (})^ ,  =0.  (/  =  0, 1 . k)\  (12 ) 

u,,  =  0.  <t  =  0.l . fe-1);  u,,  =  -^h. 

Boundary  conditions  (4),  (5),  (7)  make  it  possible  to  locate 
the  values  of  these  functions  at  the  vertical-line  nodes  i  =  0  and 
i  =  k  (cf.  Pig.  1) 


i  =  0;  u,  /  =  !  =  0-  i  =  k-,  a;,./  =  =  0; 

“*./+!  =  (2  -  ^  )  • 


(13) 


The  difference  scheme  (10),  (12),  (13)  was  used  to  perform 
on  a  computer  calculations  whose  results  are  illustrated  in  the 
diagrams. 

The  behavior  of  the  shell  was  investigated  over  the  time  of 
the  passage  of  the  longitudinal  compression  wave,  propagating  after 
the  Impact  from  the  struck  end  of  the  shell,  and  the  first  re¬ 
flected  compression  wave  spreading  from  the  fastened  end. 

The  values  of  the  parameters  entering  Into  Initial  conditions 
(3),  boundary  conditions  (7)  and  Initial  conditions  (8)  were 
taken  as:  m  -  3.64;  o*  =  0.48«10“*;  v  -  2.4;  Vo/a  -  0.0005,  0.001; 
0.002,  0.004. 

To  verify  the  computational  stability  of  the  difference 
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scheme,  the  decrease  in  the  net  mesh  h  with  time  was  calculated. 
The  initial  mesh  values  were  taken  as  >i  =  I  -  0.03.  It  turned  out 
that  a  five-fold  decrease  in  h  results  in  an  Insignificant  change 
in  the  unknown  functions  (less  than  1%) .  This  demonstrates  the 
computational  stability  of  the  difference  scheme  employed. 

To  verify  the  convergence  of  the  difference  scheme,  the  de¬ 
crease  In  the  net  meshes  h  and  I  was  calculated.  Figure  2  shows 
the  values  of  w  at  the  Instant  x  =  4.8  (reflected  compression 
wave  reaches  the  originally  struck  side  of  the  shell).  These  val¬ 
ues  were  obtained  for  Vo/a  *  0.004  with  the  following  mesh  val¬ 
ues:  h  *  I  ^  0.12  (curve  1),  h  =  I  =  0.06  (curve  2),  h  =  I  = 

=  0.03  (curve  3). 

As  can  be  seen  from  the  plot  (cf.  Pig.  2),  the  decrease  in 
net  meshes  does  not  markedly  change  the  picture  of  transverse  de¬ 
formation  of  the  shell,  which  confirms  the  convergence  of  the 
difference  scheme.  It  may  therefore  be  assumed  that  the  solution 
obtained  on  using  a  net  with  the  meshes  h  =  I  =  0.03  reflects  with 
sufficient  accuracy  the  magnitude  of  the  maximum  flexure,  the  pat¬ 
tern  of  variation  in  flexure  over  the  length  of  the  shell  and  the 
process  of  its  variation  in  time. 

The  plots  in  Figs.  3,  4,5  present  the  flexure  function  u  in 
relation  to  the  coordinate  a  at  the  time  Instants  when  the  longi¬ 
tudinal  compression  wave  propagated  as  far  as  one-half  of  the 
shell's  length  (Pig.  3,  t  -  1.2)  and  throughout  the  shell  (Pig.  4, 
T  •  2.4)  and  at  the  time  Instant  when  the  reflected  compression 
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wave  reached  the  originally  struck  edge  of  the  shell  (Fig.  5»  t  ■ 

*  ^.8). 
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Curves  1,  2,  3,  ^  correspond  to  the  following  values  of  the 
dimensionless  velocities  of  the  striking  body:  7o/a  =  0.0005, 
0.001,  0.002,  0.004. 

As  can  be  seen  from  these  plots,  during  the  passage  of  the 
first  longitudinal  compression  wave  the  greatest  growth  in  trans¬ 
verse  displacements  u  Is  observed  In  the  neighborhood  of  the  orig¬ 
inally  struck  edge  of  the  shell.  Once  the  longitudinal  wave  gets 
reflected  and  begins  to  propagate  from  the  supported  side  of  the 
shell,  however,  the  transverse  deformations  at  this  latter  side 
grow  at  a  faster  rate .  The  local  nature  of  Intense  transverse  dis¬ 
placements  Is  the  more  distinct  the  greater  the  velocity  of  the 
striking  body  Is . 

Thus,  the  qualitative  picture  of  deformation  of  the  shell 
following  a  longitudinal  Impact  Is  as  follows :  as  the  longitu¬ 
dinal  compression  wave  following  the  Impact  propagates  through¬ 
out  the  shell.  It  brings  transverse  deformations  In  Its  wake, 
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and  the  region  of  maximum  transverse  displacements  encompasses  the 
part  of  the  shell  adjoining  the  originally  struck  edge.  This  region 
becomes  Increasingly  distinct  with  Increase  In  the  dimensionless 
velocity  of  the  striking  body. 

Throughout  the  second  stage  of  passage  of  the  longitudinal 
compression  wave  transverse  buckling  decreases  In  the  neighborhood 
of  the  originally  struck  edge  and  increases  In  the  neighborhood  of 
the  supported  edge.  This  occurs  owing  to  the  superposition  of  the 
reflected  compression  wave  on  the  direct  compression  wave  at  the 
supported  end  and  hence  also  owing  to  the  attendant  marked  In¬ 
crease  In  longitudinal  stresses  at  the  supported  end.  Toward  the 
instant  when  the  reflected  compression  wave  reaches  the  original¬ 
ly  struck  edge  of  the  shell,  a  region  of  maximum  transverse  dis¬ 
placement  occurs  at  the  supported  edge,  by  analogy  with  the  pic¬ 
ture  that  prevailed  during  the  propagation  of  the  direct  compres¬ 
sion  wave. 

The  formation  of  local  regions  of  maximum  transverse  dis¬ 
placements  adjoining  either  one  of  the  sides  of  the  shell  Is  in 
qualitative  agreement  with  the  postulates  of  the  phenomenological 
theory  of  A.  Koppa  [6]  based  on  experimental  findings. 
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